Drag Reduction by a Linear Viscosity Profile 
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Drag reduction by polymers in turbulent flows raises an apparent contradiction: the stretching 
of the polymers must increase the viscosity, so why is the drag reduced? A recent theory proposed 
that drag reduction in agreement with experiments is consistent with the effective viscosity growing 
linearly with the distance from the wall. With this self consistent solution the reduction in the 
Reynolds stress overwhelms the increase in viscous drag. In this Letter we show, using Direct 
Numerical Simulations, that a linear viscosity profile indeed reduces the drag in agreement with the 
theory and in close correspondence with direct simulations of the FENE-P model at the same flow 
conditions. 



The addition of few tens of parts per million (by 
weight) of long-chain polymers to turbulent fluid flows 
in channels or pipes can bring about a reduction of the 
friction drag by up to 80% Q, II II In spite of a 
large amount of experimental and simulational data, the 
fundamental mechanism has remained under debate for 
a long time 0, 01 El ■ Since polymers tend to stretch in 
a turbulent flow, increasing thus the bulk viscosity, it 
appears contradictory that they should reduce the drag. 
There must exist a mechanism that compensates for the 
increased viscosity. Indeed, drag is caused by two rea- 
sons, one viscous, and the other inertial, related to the 
momentum flux from the bulk to the wall. For a fixed 
rate (per unit mass) of momentum generated by the pres- 
sure gradient, reducing the momentum flux can reduce 
the drag. In a recent theory of drag reduction in wall 
turbulence it was proposed that the polymer stretch- 
ing gives rise to a self-consistent effective viscosity that 
increases with the distance from the wall. Such a pro- 
file reduces the Reynolds stress (i.e. the momentum flux 
to the wall) more than it increases the viscous drag; the 
result is drag reduction. The aim of this Letter is to sub- 
stantiate this mechanism for drag reduction on the basis 
of Direct Numerical Simulations. 

The onset of turbulence in channel or pipe flows in- 
creases the drag dramatically. For Newtonian flows (in 
which the kinematic viscosity is constant) the momen- 
tum flux is dominated by the so-called Reynolds stress, 
leading to a logarithmic (von-Karman) dependence of the 
mean velocity on the distance from the wall 8]. However, 
with polymers, the drag reduction entails a change in 
the von-Karman log law such that a much higher mean 
velocity is achieved. In particular, for high concentra- 
tions of polymers, a regime of maximum drag reduction 
is attained (the "MDR asymptote"), independent of the 
chemical identity of the polymer |2|, see Fig. 2] In a 
recent theoretical paper [7j the fundamental mechanism 
for this phenomenon was elucidated: while momentum is 
produced at a fixed rate by the forcing, polymer stretch- 
ing results in a suppression of the momentum flux from 
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FIG. 1: Mean velocity profiles as a function of the distance 
from the wall [in "wall" units cf. Eq. JSJ]. The solid 
line (numerical simulations |Tl|l an d the experimental points 
(open circles) 1 1 2i represent the Newtonian results with von- 
Karman's log law of the wall observed for y + > 20. The red 
data points (squares) represent the MDR asymptote. The 
dashed red curve represents the theoretical universal MDR 
asymptote, Eq. (3) 



the bulk to the wall. Accordingly the mean velocity in 
the channel must increase. It was shown that when the 
concentration of the polymers is large enough there exists 
a new logarithmic law for the mean velocity with a slope 
that fits existing numerical and experimental data. The 
law is universal, thus explaining the MDR asymptote. 

To see how this mechanism works, consider the mod- 
ified Navier-Stokes (NS) equation for the polymer solu- 
tions : 



dU/dt + U ■ VU = -V P + V • T + is Q V 2 U 



(1) 



where vq is the kinematic viscosity of the carrier fluid 
and 7" is the extra stress tensor that is due to the poly- 
mer. Denoting the polymer end-to-end vector distance 
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(normalized by its equilibrium value) as r, the average 



dimensionless extension tensor TZ is TZi 



(nrj), and 



the extra stress tensor is (with = dUi/dxj), 

t = v p {uj - n + n- u T - dm/ dt - u ■ vn) 



(2) 



Here v p (proportional to the polymer concentration) is 
the polymeric contribution to the viscosity in the limit of 
zero shear. Note that expression (J2J for the extra stress 
tensor is valid for any dumbbell model of the polymeric 
molecule (i.e., FENE-P, Hookean, etc.). 

In 0, two simplifying approximations led to a trans- 
parent semi-quantitative theory of drag reduction. The 
first approximation is to ignore the fluctuations of 1Z as 
compared to its mean and to take 1Z sa (TV) in Eq. J3J). 
Then the main contribution of the stress tensor T in the 
modified NS equation (JJJ can be formally written in the 
form of some effective, 7?.-dependent tensorial viscosity. 
The tensorial structure of the effective viscosity describes 
different damping of different components of the velocity 
fluctuations due to the preferred orientation of polymer 
molecules along the mean shear. The second approxi- 
mation is to ignore the preferential orientation and to 



replace TZi 



TZSij . These two approximations allow 



one to simplify the dumbbell model QJ to a modified NS 
equation with an effective viscosity, vq =>• i>, and pres- 
sure, p P : 



dUi/dt + UjVjU, = -ViP + VjXwy + Uji)] , 
u = Vq + v p TI, P = p + dv/dt + U-Vv. 



(3) 



Obviously, the polymer elongation, 1Z, depends on the 
distance from the wall, leading to the corresponding de- 
pendence of the effective viscosity. Notice that the above 
approximations are uncontrolled. Their verification is 
one of the main goals of this Letter. We demonstrate 
that the dynamical model Eq. J3J) contains the essential 
properties of the full FENE-P model Eqs. 0, ©. 

In Ref. we modified the Reynolds closure approach 
to the situation in which the scalar viscosity cannot be 
neglected. This approach, that was justified by consider- 
ing a reasonable model of the coil-stretch transition of the 
polymers, resulted in the analytical form of the universal 
MDR logarithmic profile 

V+(y+) = — In (e Kv y+) , k v « 0.09 , (4) 

K v 

presented in the wall units 

TZe T = L^L/vo , y+ = yVe T /L , V+ = V/y/p 7 ! . (5) 

Here p' = —dp/dx is fixed pressure gradients in the 
streamwise direction x, L is the half- width of the channel 
(in the wall- normal direction y), and 7Ze T is the friction 
Reynolds number. Equation (2J is in excellent agreement 
with the MDR asymptote, see Fig.l. Another conclusion 
of is that in the MDR regime the polymer extension 



TZ(y) is self-adjusted in order to provide a universal linear 
profile of the effective viscosity 

^.DH^^^vV^r ( 6 ) 

Needless to say, in the viscous sublayer the viscosity is 
Newtonian: v = v$. It should be noted that the possibil- 
ity of drag reduction by increasing of the viscosity looks 
somewhat paradoxical: in the usual Newtonian case with 
constant viscosity, the drag is monotonically increasing 
with the viscosity. The point is that for the polymer so- 
lutions the effective viscosity is not constant anymore, it 
increases linearly with the distance from the wall accord- 
ing to Eq. ©. This point is the essential difference of 
the theory 7] from all previous "viscous" theories of drag 
reduction (see, e.g., |{|). Therefore a crucial test of the 
theory is to introduce such a linear viscosity profile to 
the NS Eq. J2J by hand, and see whether we observe drag 
reduction together with its various statistical aspects. 

To this aim we simulate the effective NS Eq. Q with 
proper viscosity profiles (discussed below) and show that 
the results are in semi-quantitative agreement with the 
corresponding full FENE-P DNS. All simulations were 
done in a domain 2~kL x2Lx 1.2tt L, with periodic bound- 
ary conditions in the streamwise and spanwise directions, 
and with no slip conditions on the walls that were sep- 
arated by 2L in the wall-normal direction. An imposed 
mass flux and the same Newtonian initial conditions were 
used. The Reynolds number IZe (computed with the cen- 
terline velocity) was 6000 in all the runs. The grid res- 
olution is 96 x 129 x 96 for the linear viscosity profile 
runs and 96 x 193 x 96 for the FENE-P run. The lat- 
ter was done with De T = 52.7, n p — 0.1, TZ 2 max — 1000. 
For a definition of these parameters and details of the 
numerical procedure see Ref. 0] . 

The y dependence of the scalar effective viscosity was 
close to being piece- wise linear along the channel height, 
namely v = v§ for y < y%, a linear portion with a pre- 
scribed slope for y\ < y < j/2, and again a constant value 
for 2/ 2 < V < L. For numerical stability this profile was 
smoothed out according to the differential equation 
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integrated with initial conditions ^(0) = vq, and v'{L) = 
0. We chose a = 0.04L, Vl = L/C, y 2 = 3L/4, while C 
is the dimensionless value of the slope. Examples of four 
such profiles are shown in Fig. [21 The slope of the lin- 
ear part of the viscosity profile [i.e., parameter k v in Ipjjl] 
was varied for different runs and was smaller than the 
theoretically predicted value for the MDR regime. The 
simulations with v MU n{y) require larger IZe numbers to 
sustain the turbulence. Included in the figure is the flat 
viscosity profile of the standard Newtonian flow. Since we 
keep the throughput constant, the runs differ in the val- 
ues of friction Reynolds numbers 7Ze T = ^/t^L/vq where 
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FIG. 2: The Newtonian viscosity profile and four examples 
of close to linear viscosity profiles employed in the numerical 
simulations. Solid black line: run N, solid blue: run R, red: 
run S, green: T, violet: run U. 
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FIG. 3: The reduced mean velocity as a function of the re- 
duced distance from the wall. The line types correspond to 
those used in Fig. |^| 



t w is the average friction at the wall t w = vq dU / dy. The 
decreased value of Re T is a manifestation of the drag re- 
duction measured in percentage as DR% — (t~ — t w)/ t w- 
The normalized slopes, the value of 7?e T and the percent- 
age of drag reduction for these runs are summarized in 
Table I. 

In Fig. OHwe show the resulting profiles of V^{y) vs. 
y + . The line types are chosen to correspond to those used 
in Fig. |21 The decrease of the drag with the increase of 
the slope of the viscosity profiles is obvious. Since the 
slopes of the viscosity profiles are smaller than needed to 
achieve the MDR asymptote for the corresponding 7?e r , 
the drag reduction occurs only in the near-wall region 
and the Newtonian plugs are clearly visible. 

It is most interesting to compare the effect of the lin- 
ear viscosity profile to the simulation of the FENE-P 
model in which both the throughput and TZe T are the 
same. Such a comparison was performed for the "S" vis- 
cous run, for which 7?e r = 214 and the FENE-P run 
with 7te T = 212.5. The results are presented in the 
two panels of Fig. In the upper panel the mean ve- 
locity profiles (symbols for FENE-P and dashed line for 
the linear viscosity profile) are seen to correspond very 
closely. The region with increased slope of the mean ve- 



TABLE I: DNS parameters for effective viscosity runs. 
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DR% 
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245 
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227 


13.8 


0.035 
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214 


21.6 


0.042 
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10 


197 


36.9 


0.051 
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12 


185 


42.0 


0.065 





0.4 y/L 0.6 



FIG. 4: Upper panel: the reduced mean velocity as a function 
of the reduced distance from the wall. Lower panel: Reynolds 
stresses across the channel. Continuous line: Newtonian. 
Dashed line: linear viscosity profile. Symbols: FENE-P. 
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locity fluctuations. The increase in r.m.s. streamwise 
(Vj~) and decrease in r.m.s. wall-normal (Vj~) velocity 
fluctuations were observed in many experiments and nu- 
merical simulations of drag-reducing flows. In Fig. [5] 
we compare these quantities for FENE-P and "S" runs. 
Clearly, the correct trend of the r.m.s velocity fluctua- 
tions is observed, indicating that the important features 
of the mechanism of drag reduction are correctly cap- 
tured by the model. An almost quantitative agreement 
is reached in the region where drag reduction actually 
occurs ( y/L w 0.1 — 0.3). 

In conclusion, we showed that the simple linear viscos- 
ity model ©, © faithfully demonstrates drag-reducing 
properties and, surprisingly enough, the amount of drag 
reduction increases with the increase of the slope of the 
viscosity profile. Even more interestingly, the behavior 
of objects like the Reynolds stress or the velocity fluctu- 
ations in the elastic sublayer are in close correspondence 
with the full FENE-P model, indicating that the mecha- 
nisms of drag reduction proposed in operates similarly 
in both cases. 

This work was supported in part by the US-Israel BSF, 
The ISF administered by the Israeli Academy of Science, 
the European Commission under a TMR grant and the 
Minerva Foundation, Munich, Germany. 
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FIG. 5: The rms streamwise ( upper panel) and wall-normal 
(lower panel) velocity fluctuations across the channel. The 
line types correspond to those used in Fig. 0] 



locity (e.g., the region where the drag reduction occurs) 
is limited to approximately y + < 50, or, in natural units, 
to y/L < 0.25. In the outer region (Newtonian plugs, 
y + > 50 or y/L > 0.25) the polymer molecules in the 
FENE-P model are supposedly unstretched and the poly- 
meric contribution to the effective viscosity is small. In 
contrast, in the viscous model the effective viscosity is 
maximum in this region and may be much larger than 
vq. Therefore, if the mechanism of the drag reduction 
in viscous model is the same as in full FENE-P model, 
we expect that all statistical quantities qualitatively co- 
incide for both models in the elastic region y/L < 0.25, 
but may differ in the Newtonian plugs y/L > 0.25. 

In Fig. 01 lower panel, we show the normalized 
Reynolds stresses. Clearly, in the elastic region y/L < 
0.25 both drag-reducing models coincide; this is a strong 
evidence that the reduced model @ captures all essen- 
tial properties of the full model Q and the mechanisms 
of drag reduction are the same in both cases. 

Another important characteristic of drag-reducing 
flows is behavior of the root-mean-square (r.m.s.) ve- 
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